Dirac spinors in solenoidal field and self adjoint extensions of its Hamiltonian 
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We discuss Dirac equation (DE) and its solution in presence of solenoid (infinitely long) field in 
(3+1) dimensions. Starting with a very restricted domain for the Hamiltonian, we show that a 
1-parameter family of self adjoint extensions (SAE) are necessary to make sure the correct evolution 
of the Dirac spinors. Within the extended domain bound state (BS) and scattering state (SS) 
solutions are obtained. We argue that the existence of bound state in such system is basically due 
the breaking of classical scaling symmetry by the quantization procedure. A remarkable effect of 
the scaling anomaly is that it puts an open bound on both sides of the Dirac sea, i.e., E 6 (-M, M) 
for v 2 [0, 1)! We also study the issue of relationship between scattering state and bound state in the 
region v 2 £ [0, 1) and recovered the bound state solution and eigenvalue from the scattering state 
solution. 

PACS numbers: 03.65.-w; 02.30.Sa; 03.65.Pm 



I. INTRODUCTION 

The problem of a Dirac particle in magnetic fields 
background has been studied extensively in literature 
P, H, H, 0] • The solution in a uniform field background [l| 
is known for a long time. For a solenoidal background, 
the Dirac equation has been solved in (2+1) dimensions 
@, B 0|- -"-t shows that a self-adjoint extension of the 
Dirac Hamiltonian is necessary to obtain the solutions. 
In (3+1) dimensions, self-adjoint extension of the Dirac 
Hamiltonian has also been studied but it involves a 
(5-sphere potential rather than solenoid interaction. In 
Ref. Q, self-adjoint extension of Dirac Hamiltonian in 
solenoid field and an uniform magnetic field B has be 
studied in both (2+1) and (3+1) dimensions elegantly. 

The consequences of considering self-adjoint extensions 
for the problem of Dirac particle in solenoidal field is 
recog nized in Ref. Q. Because SAE @, ©, G3 is a 
method which allows to consistently build the all pos- 
sible boundary conditions under which the Hamiltonian 
is self-adjoint. It is therefore expected that [2| the results 
would be function of the boundary conditions, character- 
ized by parameters. 

In this paper, we study the solution of the Dirac equa- 
tion in only infinitely long solenoid field in (3 + 1) dimen- 
sions. Because, the magnetic field of an infinitely long 
solenoid is immensely important for various reasons. For 
example, it is used to study the Aharanov-Bohm [llj 
effect. Similarity with this field also helps one to un- 
derstand the physics of cosmic strings [12] ■ Although, in 
principle problem in Ref. Q should reduce to the problem 
involving only solenoid field by making the uniform mag- 
netic field B — 0. But it is a nontrivial task. So it is im- 
portant to study Dirac free Hamiltonian in only solenoid 
field background in (3+1) dimensions and solve it with 
the most generalized boundary conditions, such that the 
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Hamiltonian is self-adjoint. Some comments about the 
difference between method of [6] and ours will be made 
in Sec. VII. 

In order to get the most generalized boundary con- 
dition for our problem, we start with a very restricted 
domain. Obviously the Hamiltonian is not self-adjoint 
in that domain. We then go for a self-adjoint exten- 
sions of the Hamiltonian by using the von Neumann's 
method [l5(. We map the problem from the Hilbert 
space Li(drd(\)dz;C^) to L2(drd<f>dz;C) keeping in mind 
that we seek self-adjointness of the whole (4 x 4) Dirac 
Hamiltonian in solenoid field background. In the Hilbert 
space Li(drd$dz; C) we get a well known one dimensional 
Schrodinger eigenvalue problem of inverse square interac- 
tion. Inverse square interaction is known to be classically 
scale invariant. But, quantization of this system shows 
that the classical scale symmetry is broken for some val- 
ues of the self-adjoint extension parameter and for a very 
limited range —1/4 < g < 3/4 of the coupling constant g 
of the inverse square interaction. Breaking of this scaling 
symmetry by quantization is known as scaling anomaly. 
The indication of scaling anomaly in the system is man- 
ifested through the formation of bound state, which is 
supposed to be absent if scaling symmetry is present even 
after quantization. In case of Dirac particle in solenoid 
field, this scaling anomaly is also responsible for the for- 
mation of bound state. It is also possible to comment on 
the bounds (upper and lower bounds) of the bound state 
energy of the Dirac particle on solenoid filed background. 
By taking appropriate limit, it is even possible to find the 
upper and lower bounds of the energy of the free Dirac 
particle. 

The paper is organized as follows: In Sec. II, we 
state the problem of Dirac equation (DE) in solenoid field 
background. In Sec. Ill, we state the usual symmetric 
boundary condition using regularity and square integra- 
bility argument and we perform the required self adjoint 
extension. In Sec. IV, we discuss the bound states of the 
radial Hamiltonian and we get the bound state condition. 
In Sec. V, we discuss the scattering states of the radial 
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Hamiltonian and get the corresponding condition for it. 
In Sec. VI, we discuss the scaling anomaly present in the 
system. We discuss in Sec. VII. 



II. SOLUTIONS OF THE DE IN SOLENOID 
FIELD 

We consider a Dirac particle of mass M and charge eQ 
in solenoid field background. Due to the cylindrical sym- 
metry of the problem it is better to use the cylindrical 
polar coordinate system, where the spatial coordinates 
are denoted by r, <p, z. In any orthogonal coordinate sys- 
tem, the Dirac equation for a particle with charge eQ and 
mass M can be written as Ref. 



[y (iD„ - eQA, 



M 



-Di log (hi/12/13/fti) 



0, 



(1) 



where h^s are scale factors of the corresponding coor- 
dinate satem. In our case since we are using cylindrical 
coordinate system, h% = 1, /12 = r, /13 = 1. the derivative 
Dfj, = (h fJl )~ 1 d fl , where no summation over \i is implied 
here and is the infinitely long solenoid vector poten- 
tial, A = —<b. One can make a conformal transformation 

r 1' r T 

|13l |. which reduces ([1]) to simpler form with zero spin 
connection as 1131 . 



Da - eQA^ 



M 



*(t,r,0,*) = O, 



(2) 



where the relation between 'Jv and W is \& r = 
(l/\/r) exp( — gi ^Ss)^. S3 is the generator for rotation 
along z-axis |13| | . The normalization condition of the new 
wave-function "J is 



drd(f)dz r, 0, r, 0, z) 



1. 



(3) 



To solve the Dirac equation §2§ we take the trial solution 
of the Dirac equation of the form 

*(t, r, 0, z) = e -«« e -<m* e -ip.» ( ) , (4) 

where and x are 2-component objects. We use the 
Pauli-Dirac representation of the Dirac matrices 



7 



a 1 



-a 1 



7 = 



1 

-1 



(5) 



where each block represents a 2 x 2 matrix, and <7jS are 
the Pauli matrices. Note that the same 7^ matrices are 
used which we use in Cartesian co-ordinate system. For 
detail discussion about the Dirac equation in cylindrical 
co-ordinates and the form of the 7^ matrices for cylin- 
drical co-ordinate system see Ref. [l3|]. Multiplying Eq. 
@ by 7 from left, and using Eqs. (U]) and |(SJ| we get 



E — M A 
A E + M 



0(r) 
X(r) 



= 0, 



(0) 



where A = cr. (iD^ — eQA^j = i^d, 
Eq. ^ is divided into two coupled equations 



^-^cr 2 +a 3 p z 



x(r) 



ia x d r 



m—eQa ^2 



a 3 p z 



E-M 



ia x d r 



a 3 p z 



E + M 
Eliminating x( r )> we obtain 



m — eQa ^2 



a 3 p z 



E 2 - AP 



(7) 
(8) 

(9) 



There will be two independent solutions for <f>(r), which 
can be taken, without any loss of generality, to be the 
eigen-states of a z with eigenvalues s = ±1. This means 
that we can choose two independent solutions of the form 



F+(r) 






F_(r) 



Since a 3 



r ), Eq. ([9]) becomes 



") 1.= 



±i = 



dr 2 



(m-eQa^l/2) -1/4 



E 2 - M 2 



(10) 



") (11) 



For s = + 1, using (fT0|) in (fTTj) . the differential equation 
satisfied by F + is 



d 2 F + 



+ A - 



v 2 - 1/4 



F+ = 0, 



(12) 



where 

A = (E 2 — M 2 - p 2 z ) 1/2 , v = m- eQa - 1/2 . (13) 

It should be noted that, (fl"2"|) can be considered as a well 
known one dimensional Schrodinger eigen- value equation 
with inverse square interaction yl, 0, HO] ■ It shows clas- 
sical scale symmetry and scaling anomaly, which will be 
discussed separately in Sec. VI. Now we seek the solution 
of this equation, which is of the form F + = r 2C to v {\r). 
Where denotes J, Y, H 1 , H 2 or any linear combination 
of these functions with constant coefficients. From |8]), 
we find the two lower components and write the spinor 
as 



U+(r,p z ) 







Vz 



E + M 
iX 



\ E + M 



(14) 
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where the normalization has not been specified. For s — 
— 1, we get the same Bessel differential equation (fT2"j). 
with v replaced by v + 1. The solution, which can be 
obtained similarly for this case also, is of the form 



U-(r,p M ) = 





- iX 



\ 



E + M 
P 



(15) 



A similar procedure can be adopted for negative fre- 
quency spinors. In this case, it is easier to start with 
the two lower components first and then find the upper 
components. The negative energy spinors are found to 
be 



V+(r,p z ) 



Pz 



E + M 
iX 



r»%(Ar) 



E + M 

\ o 

/ -iX 



rW v+l {Xr) 



(16) 



V-(r, Pz ) 



E + M 



r*Sf„(Ar) 



E + M 



r^ v+l {Xr) 







(17) 



V rW v+1 {Xr) J 



III. SAE OF RADIAL HAMILTONIAN 



However, in domain D{H) the radial Hamiltonian H(r) is 
not self-adjoint. A symmetric Hamiltonian is self-adjoint 
if its domain coincides with that of the domain of its ad- 
joint, i.e, D(H) = D(W). The condition at the origin 
makes the Hamiltonian non self-adjoint. The domain of 
the adjoint Hamiltonian H^(r) is D(H^) — ip(r), where 
4>{r) G L 2 (dr; C 4 ). We see that D(H) ^ D(H^), indicat- 
ing that H (r) is not self-adjoint. To make the Hamilto- 
nian self adjoint [l5| we use von Neumann's method of 
deficiency indices. It requires the construction of eigen- 
space of with eigenvalue +iM (M ^ is inserted 
for dimensional reason). The up spinors (particle state) 
for the eigen-space Z? ± are 



= N 



r2Hy(e ±l iXir) 







M(fV^' 2(e±,;iAir) 



(21) 



where Ai = (2M 2 +p^j 2 and N is the normalization 
constant. Here we have written cjx^ for up spinor (particle 
state) only. Similarly we will get (j> for all other spinors. 
Since we will do our calculations for spin up particle state 
only, (|2"T|) is sufficient for us. Looking at the asymptotic 
form of the Hankcl functions 



Hl{z) 



Hl(z) 



2 

7TZ 



1/2 



for (—Ti < argz < 2tt) 



2 

TTZ 



1/2 



for (— 2ir < argz < it) 



(22) 



(23) 



We now discuss the self-adjointness problem of our 
system. Therefore we need to know the radial eigen- 
value problem, which can be obtained from ([2|). It can 
be shown that the radial eigenvalue equation is of the 
following form 



H(r)S(r) = ES(r) , 



(18) 



where the radial Hamiltonian and the eigenfunction are 
given by 



H(r) 



M -A 
-A -M 



.5(r) 



<P\r) 

x( r ) 



(19) 



where A = ia x d r + m ~ eQa cr 2 + a 3 p z . The differential 
operator H(r) is symmetric over the domain D(H) — 
ip(r), where ip(r) £ L^idr; C 4 ) and ip(0) = 0. This means 
that for ipi,ip2 € D(H), the radial Hamiltonian H(r) 
satisfies the condition 

/>OC />OC 

/ dril){(r)H(r)ip 2 (r) = / dr [ J ff(r)V'i(r-)] t Va(r)(20) 
Jo Jo 



we find that the upper end of the integrals for evaluating 
the norms of ^ are finite for any v. However, near r = 0, 
the Hankel function behavior can be found from the short 
distance behavior of the Bessel function 



J v {z) 



2T(l + i/)' 



(1/54-1,-2,-3, 



(24) 



Considering all components of the spinor of (|2"Tj) , we find 
that (j)^ are square integrable only in the interval 



< v l < 1. 



(25) 



Since beyond this range there is no square integrable so- 
lutions <ft , the deficiency indices, which are dimensions 
of the eigen-space D ± . 



are zero , i.e, 



n± = dim(D 1 ), 



n + 



(26) 



(27) 
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The closure of the operator H (r) is the self-adjoint exten- 
sion for the case (|2"T)) . We therefore concentrate for the 
interval (|25| to carry out self-adjoint extensions. The 
existence of complex eigenvalues for H'(r) emphasizes 
the lack of self-adjointness. The self-adjoint extensions 
of H(r) are labeled by the isometries D + — ► D~ , which 
can be parameterized by 



0+(r) 



(28) 



The correct domain ( It is better to call projection of the 
domain on spin up particle state direction rather than 
only domain because we are concentrating on spin up 
particle state only here ) for the self-adjoint extension 
H u (r) of H(r) is then given by 



D U (H") ee D(H) + <P + (r) + e ia >-(r) , 



(29) 



where cu £ R(mod27r). In the next section we find out 
the bound state solutions using the domain D UJ (H UJ ). 



IV. SOLUTIONS OF RADIAL HAMILTONIAN 

From the trial solution (0| it is clear that the spinors 
along the z-direction is free, as it should be, because there 
is no constraint in the z-direction. We therefore try to 
investigate whether the spinors are bound in the radial 
direction. From now on bound state means bound in the 
radial direction. Throughout our calculation, we will use 
the spin up state. Calculation for all other spinor states 
are similar. From the general spin up state (114|) it is easy 
to see that it serves as a square integrable spin up state if 
weuse^(Ar) = Z#(Ar) and if A = {E 2 -M 2 -p 2 z )i = iq, 
where q is real positive. Similarly we could have used 
^(Ar) = H 2 (Xr) as square integrable function if we took 
A = (E 2 — M 2 — p 2 ) 2 = —iq, where q is real positive. We 
will use Hankel function of the first kind ff* to express 
bound state solution. The bound state spinor of spin up 
particle is then found to be 



/ r*#i(Ar) \ 




U+(r,p z ) = BVE + M 



Pz 



E + M 
iX 



V E + M 



r^Hl +1 {Xr) 



(30) 



where B is the normalization constant. To find out bound 
state eigenvalue we have to match the limiting value r — > 
of spinor (|30[) with the limiting value r — > of the 
domain (f2"9")l . For the sake of simplicity, we set p z = 
before matching at the origin. The relevant range of v 
is given in (|25p . In this range, the leading r-behavior of 



the domain (|29| for small r is given by 

/ A(A!)^ +1 /2 \ 





ijj(r) +4> + {r) +e^<p-(r) 



where ^4(Ai) 



N 



-Ni^- 



sin u-k 2"r(l+jy) 





V D(A x )r- i/ - 1 / 2 
K_ 



> (31) 



D(Ai) _ --tvt ain( „ +1)7r 2 -,-i r( _^) 
and the leading term of the spinor (|3"0")) is 



U+(r,p z ) 



( A{X)r»+ 1 / 2 \ 



V £>(A)r-^ 1 /2 J 

K. p -tvi fj(X) 



(32) 



M- 



where A(X) = B^E , ^ sin _ 2 „ 

-BiVE^M ^±^ ■ Since u <r>P>) e 

Z? w (iJ"), the coefficients of leading powers of r in (j3"T|) 
and ([32]) must match and this gives the eigenvalue equa- 
tion 



_2^+i/2_ 



sin (lo/2 + ttv/2) 
sin (ui/2 + nv/2 + tt/4) 



(33) 



The left hand side of (|3"3"|) is positive. So right hand 
side should be positive and to ensure that we impose the 
condition cot(^ + ^) < —1, which is the bound state 
condition. Similarly we may get all other spinors for 
the bound states. We move to the next section for the 
discussion of scattering state (SS) solutions. 



V. SS SOLUTIONS OF RADIAL HAMILTONIAN 

The scattering state spinor of spin up particle is 

/ r\A v \ 




U+(r, Pz ) = BVE + M 



Pz 



E + M 
iX 



(34) 



V ETM r ~ 2A » +1 J 

where B is the normalization constant, 
A v = [a{X)J v {Xr)+b{X)J- v {\r)] and A u+1 = 

a(A)J i/+ i(Ar)-6(A)J_ 1 ,_ 1 (Ar)], a(A), 6(A), 5(A) 

and 6(A) are constant coefficients. To find out eigen- 
value for the scattering state we have to match the 
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limiting value r — > of the spinor ([34)1 with ([29)) . For 
simplicity of calculation we set p z = in (|34[) and 
(|29p before matching. In the limit r — > 0, the spinor 
([3"4")l looks like (|32[) but now the coefficients of differ- 
ent powers of r are, A(X) = B\/E + Ma(X)- 

x—- 1 



that, we rewrite (|12[) in the form of time independent 
Schrodinger equation 



A" 



and the limit r — > 



D(A) = -Bzv^TM6(A) 2 _,_ lr( _ y) 
of is given in (|3T|) . Again equating the respective 
coefficients and comparing between them we get the 
eigenvalue equation 

(E + M) 1/2 (A) ( X_^ 2u+l 
Ai , 



(25 



M) 1/2 6(A) 

sin (w/2 + 7w/2) 
sin (uj/2 + ttv/2 + tt/4) 



(35) 



The right hand side of (|35[) has to be positive in order to 
get scattering state solutions. Similarly we may get all 
other spinors for scattering states. We now want to show 
the relation between scattering state and bound state 
[1 , [TtJ for completeness of our calculation. To calculate 
that we expand (|34[) in the limit r — > oo. The leading 
term in the asymptotic expansion of U+(r,p z ), without 
normalization is given by 



U+(r, Pz ) 



A(X)e iXr + B(X)e 




— i\r 



\ 



+ ETM B ^ XT 
C{X)e lXr + D(X)e~ lXr 



(36) 



where 
A(A) 



different 



E+M 
2ir\ 



E+M 
2tt\ 



E — M 
2ir\ 



E-M 
2ir\ 



coefficients are found to be 

a(A)e- l f ^ +1 /2) + & ( A ) e if (^-1/2)] ; 

a (A) e 'f (-+V2) + 6 ( A ) e -if (^-i/2)j ^ 

a ( A ) e -if(^+3/2) _ 5( A ) e *f (u+1/2) 

6(A)e^+ 3 / 2 ) -b{\)e-^^ +1 /V 



B(X) = 
C(A) = 
D(X) = 

Now it is easy to see from the asymptotic expansion ([36)) 
that e~ lXr blows up on the positive imaginary A-axis 
but the other part e lXr decays exponentially there. So 
it reasonable to set the coefficients B(X) and -D(A) zero 
for purely positive imaginary A to get square integrable 
behavior at oo. This corresponds to bound state. Using 
([33]) and ([56]) it can be shown that the bound state 
eigenvalue is again given by (|33[) . Similarly it can also 
be shown that for purely negative imaginary A it is 
reasonable to set the coefficients A(X) and C(A) zero in 
order to get bound state solution. 



VI. IMPLICATION OF SCALING ANOMALY 

As pointed out in section II, we now discuss the clas- 
sical scale symmetry and scaling anomaly of (|12[) . To do 



H r F s (r)=£F s (r) 



(37) 



where the Hamiltonian H r — ~d 2 — {v 2 — 1/4) /r 2 , v = 
m — eQa — 1/2 and the eigen- value £ = E 2 — M 2 . Due 
to the cylindrical symmetry we just consider the problem 
on x-y plane, by setting p z = 0. This is a well known one 
dimensional inverse square problem, discussed in many 
areas of physics, from molecular physics to black hole [la, 

El, IlllmMl^lllllll^lillllliilii- Classically 
(13"T1) is scale covariant. The scale transformation r — > 
£r and t — > £ 2 t (£ is scaling factor) transforms the the 
Hamiltonian Ti r — -> (l/£ 2 )7i r . Scale covariance of 7i r 
implies that it should not have any bound state. But, 
it is well known that there exists a 1-parameter family 
self-adjoint extensions (SAE) of H r and due to this SAE, 
a single bound state is formed. The bound state energy 
for v 2 € [0, 1) is given by 



E 2 - M 



'sin(E/2 + 37rz//4) 
sin(E/2 + 7r^/2) 



(38) 



where E is the self- adjoint extension parameter for TL r . 
Existence of this bound state immediately breaks the 
scale symmetry, which leads to scaling anomaly. 

If we assume that the energy of the Dirac particle E 
to be real, then from the bound state condition £ < 0, 
we get a bound on the energy of the Dirac particle to be 
E S (— M, M). This bound still holds for the free Dirac 
particle, which we get by taking limit a — > in (|38|) . 

Quantum mechanically, scale transformation is asso- 
ciated with a generator, called scaling operator A = 
\{j"Pr + Prf), where p r = —i-j-- It can be shown that 
the action of the operator A on a generic element of 
the self adjoint domain D^(Tt r ) does throw the element 
outside the domain for some values of the self-adjoint 
extension parameter E. This indicates that there is 
scaling anomaly, occurred due to self-adjoint extensions. 
However, it can be shown that for E = —vk/2 and 
E = — 3^7r/2, the action of the operator A on any ele- 
ment of the domain D^(Tl r ) does not throw it outside 
the domain. So, in these two cases scaling symmetry is 
still restored even after self-adjoint extensions and thus 
the bound state does not occur. 



VII. CONCLUSION AND DISCUSSION 

In this paper we calculated the solution of the Dirac 
equation in the field of an infinitely long solenoid. We 
showed that there is nontrivial bound state and as well 
as scattering state solution in the range v 2 g [0, 1). We 
showed only spin up particle state details of self-adjoint 
extensions in our calculation, but for all other spinors 
calculations are similar. We point out that for v = 0, the 
solution of Wcfr 1 = zLiMcj)^ involves r^HKiXir), which 
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is not square integrable at the origin. So at v = de- 
ficiency indices are zero, i.e, n + = ri_ = 0, that means 
closure of domain is the self-adjoint extension for v = 0. 
We have projected the whole problem from 4-component 
Hilbert space (L2(drd<fidz;C 4: )) to 1-component Hilbcrt 
space (L2(drd<fidz;C)) keeping in mind that we make 
the whole (4 x 4) Dirac Hamiltonian self-adjoint unlike 
Ref. [6|, where they discussed the self- adjoint issue pro- 
jecting the problem on Li{rdrd$dz;C 2 ). The significant 
difference between these two approaches is that in our 
case we have a 1-parameter family of self-adjoint exten- 
sions unlike Ref. [6j] where they have 2-parameter family 
of self-adjoint extensions. As in Ref. (15 , [l7| we stud- 
ied that there is a relation between scattering state and 
bound state and it occurs for the purely imaginary value 
of A. This has been done by demanding square integrabil- 



ity of the scattering state at spatial infinity. We showed 
that scattering state eigenvalue reduces to bound state 
eigenvalue for purely imaginary A. Finally we discussed 
the implications of scaling anomaly on Dirac particle in 
background solenoid field. The energy of the Dirac par- 
ticle in solenoid field background has bounds from both 
sides E £ (— M, M). By taking appropriate limit, these 
bounds were also shown to hold for free Dirac particles. 
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